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Abstract. We consider nonlocal modification of the Einstein theory of gravity in framework of the pseudo- 
Riemannian geometry. The nonlocal term has the form where “H and are differentiable 

oo 

functions of the scalar curvature R, and !F(n) = ^ /„□" is an analytic fimction of the d'Alambert operator 

fi=0 

□. Using calculus of variations of the action functional, we derived the corresponding equations of motion. 
The variation of action is induced by variation of the gravitational field, which is the metric tensor 
Cosmological solutions are found for the case when the Ricci scalar R is constant. 


1. Introduction 

Although very successful, Einstein theory of gravity is not a final theory. There are many its modifica¬ 
tions, which are motivated by quantum gravity, string theory, astrophysics and cosmology (for a review, 
see ii)- One of very promising directions of research is nonlocal modified gravity and its applications to 
cosmology (as a review, see j^rQ] |0])- To solve cosmological Big Bang singularity, nonlocal gravity 
with replacement R ^ R + CRlF(n)R in the Einstein-Hilbert action was proposed in |@|. This nonlocal 
model is further elaborated in the series of papers SQ. 

In (l^ we introduced a new approach to nonlocal gravity given by the action 

s = jd^xyrg(j^+R-^r{a)R), ( 1 ) 

where the d'Alembert operator is □ = ff = det{gi^,v). The nonlocal term R“^;F(n)R = 

/o + R~^ /nO"R contains /o which can be cormected with the cosmological constant as /o = This 

term is also invariant under transformation R —> CR, where C is a constant, i.e. this nonlocality does not 
depend on magnitude of the scalar curvature R 0. 
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In this paper we consider n-dimensional pseudo-Riemannian manifold M with metric of signature 
(n_, n+). Our nonlocal gravity model here is larger than |[T]| and given by the action 

S = £ + 9{(R)r(D)^(R)) ^d”x, (2) 

which is a functional of metric (gravitational field) where and ^ are differentiable functions of the 
scalar curvature R, and A is cosmological constant. 


2. Variation of the action functional 

Let us introduce the following auxiliary functionals 

So= f (R-2A)^d"x, Si= f 9{(R)r(DmR)^d"x. (3) 

Jm Jm 

Then the variations of So and Si can be considered separately and the variation of ([2j can be expressed as 

6S = + 6Si. (4) 

IottG 

Note that variations of the metric tensor elements and their first derivatives are zero on the boimdary 
of manifold M, i.e. = 0, ddAyfivhm = 0. 

Lemma 2.1. Let Mbe a pseudo-Riemannian manifold. Then the following basic relations hold: 


y ‘^aa i/ ^ era' 

r|;. = ^inVW. 


Vi^i 


= -gg^iv5g>^\ 
5 a/m = £g^iv a/Mi5/^ 


5R = R, 


Lemma 2.2. 


d”x = 0. 


On the manifold M holds I g^''5Rfiv -y/M 

Jm 

Proof. Let W^' = -g^“^bTJia + g^''^^‘j!ia- Then it follows 

dx- ^ dx- 

LFsing Lemma [2.II we get 






( 5 ) 


^q^“^ dT^a dqi^'’ <?r“ 

dx^ ^ dx^ dx'' 1-“^ ^ dx'' ^ /ja V fia) v/S 


( 6 ) 














Ivan Dimitrijevic et ah/Filomat xx (yyyy), zzz-zzz 


3 


Moreover, using again Lemma|2!T]we obtain 
1 d 


^IW\ 


dx'' 


dT" 


fia 

rP 


fia 

rP 


dx'’ 


- ni^vt _ 6r“ + r'* sr" - P I 


-/M' 


^r« ^r" 

^ ^ P r« _ r« rM - 

(Jx" f^X'’ 


■r^a/?v-r;,ry = /’'5V- 


1 d 


Finally, we have 

' 'Viii*' 

Using the Gauss-Stokes theorem one obtains 

r /^(^w'')d”x= r irdu,. 


(^IT) and r /''6V\^d”x= f A(^vr)d«x. 
Jm Jm 


( 7 ) 


( 8 ) 


( 9 ) 


Since 6^^,v = 0 and bd^gf^v = 0 at the boundary dM we have W\su ~ 0. Then we have VPdcTv = 0, that 
completes the proof. □ 


Lemma 2.3. The variation of So is 

6So = r Gf,v ^J\g\6g'^''' d"x + A f d"x, 

Jm Jm 

w/iere Gf,v = Kjiv - \^0iiv is the Einstein tensor. 

Proof. The variation of So can be found as follows 

5So = f 6((R - 2A) VM) d”x = f 6(R ^/\g\) dF’x-lA f 6^/^g| d"x 
Jm Jm Jm 

= f ( + R6 x/\g\ + Ag^^y x/\g\^0^'') 

Jm 

= - ^R ^/^g|g^lvbg>^^ + Ag^,y ^/^g|bg^‘^j d"x 

= r Gf,v VMS/''d"x + A r VM6/''d^x + r /^6Rf,vVM 

Jm Jm Jm 

Using Lemma [2.21 from the last equation we obtain the variation of Sq. □ 
Lemma 2.4. For any scalar function h we have 

r hbR ^J\g\ d"x = f (hR^^y + g^yOh - V^iVyl^ bg^'^' ^J\g\ d"x. 

Jm Jm ' 

Proof. Using Lemma f2.II for any scalar function h we have 

r hbR x/\g\dPx= f (/xVb/'' + hg^.ynbg^'' - ^ d”x. 

Jm Jm ^ 


( 10 ) 


d"x. 


( 11 ) 


( 12 ) 


( 13 ) 
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The second and third term in this formula can be transformed in the following way: 


r VM d"x = r VM d"x, 

Jm Jm 

f W^V^6/^^d”x= r V^,V^h6yf^''^/iild"x. 

Jm Jm 

To prove the first of these equations we use the Stokes theorem and obtain 

r hy,,o6f-^/^i^d"x= f %MvVaV“6/"^^d”x = - f ^ d”x 

Jm Jm Jm 

= f <f^vV“V«;i6/"^d"x= r 6yf‘''x/fil d^x. 

Jm Jm 

Here we have used ^yfjf/iv - 0 and = □ fo obfain the last integral. 

To obtain the second equation we first introduce vector 


(14) 

(15) 


(16) 


N'' = hVydf'' - V^h5y^‘\ 


(17) 


From the above expression follows 


Vf.N'' = V^(W^,b/^ - WMgn = bg^'' - VyhVybg'^'' 

= hV^V,bgf‘^ - y^y,h bg^‘\ 


(18) 


Infegrating yields J^VyN^' d”x = f^^N^^n^ddM, where is the unit normal vector. Since 

^^IdM - 0 we have that the last integral is zero, which completes the proof. □ 


Lemma 2.5. Let 6 and \jj be scalar functions such that bflsM - 0- Then one has 


r dbnip ^J\g\ d"x = ^ r g^^dad dpipgftvbg^" ^J\g\d"x- f dvf6gf‘''-\/igl d"x 
Jm 2 Jjyj 

+ r aG 6f yj\g\ d"x + ^ f g^vGDfbg^'' yj\g\ d”x. 

Jm 2 Jjyi 

Proof. Since 6 and f are scalar functions such that b\p\dM = 0 we have 


r Gbaf xjigl d"x - f ddab( y]\g\f^d^f) d"x + f 06 [ J= J«( yl\g\g“^dpf) y]\g\ d"j 
Jm Jm Jm y\g\) 

= f ^a(Jb(^/^g^g‘^l^^llf)) d^x - f JaJ b(^/^g^g‘^f^Jf!f) d”x + ^ f GgyvDfbgf^''^/^g^ d”x. 
Jm Jm ^ Jm 


(19) 


( 20 ) 
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It is easy to see that da{db{ d"x = 0. From this result it follows 

I dbnxp ^J\g\ d"x - - f g“^dad dpip6{ ^/\g\) d"x - f dad dpxpdg"^ ^/\g\ d"x 

Jm Jm Jm 

- f ^l\g\dad dpbTp d"x + ^ r dgfiyOipbgf^'’ ^J\g\ d"x 

= l f VM d"x - r (9^0 dv 1/^6/''VM d”x 

- f dfiig"^ ^/^g^dad bip) d"x + f dfj(g"^ ^/^g^da9) btl> d"x 

Jm Jm 

+ l f giivdlFi^bg^''' VM d"x 

= l f g'^^dad dfi^gfivbg^"' VMd"x- f Vm d”x 

^ Jm Jm 

+ I 00 6i/> -^/i^ d"x + ^ r g^ivdnxpdg^'^ ^/\g\ d"j. 

Jm 2 Jm 


At the end we have that 


r 06ni/^ VM d"x = i r g^l^dad dp^g^vbg^'’ VM d"x 
Jm ^ Jm 

- I d^d dyi^dg^"^ ^/\g\ d"x + j ad dip ^j\g\ d”x + x ( g^vdaipbg^'' ^J\g\ d”x. 
Jm Jm 2 Jm 


□ 


Now, after this preliminary work we can get the variation of Si. 
Lemma 2.6. The variation of Si is 

^Si = -l: f J,v'K(l?)^(□)@(W^^d”x+ r (v®- 

^ Jm Jm 

_ oo «-l _ 

+ 9 E E (j,v (j“n''K(R)0«n"-i-'@(R) + □'<K(R)n"-'@(R)) 

^ n=l ;=0 dw 

- 20^n'<K(R)0va"“^“'@(R))6/^ d”x. 


( 21 ) 


( 22 ) 


if/iere = Vj,Vv - ® = ‘?^'(R)J^(n)@(R) + &'(R)T'{o)9'{{R) and' denotes derivative with respect to R. 

Proof. The variation of Si can be expressed as 

6Si = £ (<N(R)r(a)^(R)6( y^) + d(<N(R))r(n)^(R) y^ + <K(R)b(!r(n)@(R)) y^ d"x. (23) 


For fhe first two integrals in the last equation we have 

Ji= r <N(R)r(amR)^(^/ii^)d’’x = -l f (7i,v'H(R)r(n)@(R)6/"y^d”x, 

l2= f d('N(R))r(a)^(R)x/lgld"x= f <K'(R)6R r(n)@(R) y^ d”x. 

Jm Jm 


( 24 ) 
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( 25 ) 


Substituting h = 'K'(K) !F(n)@(R) in equation lfT 2 ll we obtain 

Jm 

The third integral can be presented using linear combination of the following integrals 

/„= r 9{(R)^(^"^(R))^/\i\d"x. (26) 

Jm 

Jo is the integral of the same form as h so 

Jo= f (V^'(R)'H(R)-fC,,(@'(R)<H(R)))6/’'^^d"x. ( 27 ) 

Jm 

For n > 0, we can find }„ using ifT^ . In the first step we take 6 = ^{R) and ip = □”“^^(R) and obtain 

U = \J /''<9„<K(R) ^^□”-i@(R)v 6 /" Cl”x - r <9,<W(R)5,n”-i@(R)6/"^^d”x 
^ Jm Jm 

+ r □<H(R)6n”-i@(R)^d"x + t r ^^„-H(R)□”@(R)6/''^^d"x. 

Jm ^ Jm 


(28) 


In the second step we take 6 = □'7F(R) and i/> = □” ^@(R) and get the third integral in this formula, etc. 
Using 1 IT 9 I 1 n times one obtains 


n-1 


/„ = q y r U,vd"n'<H(R)d„n"-i-'@(R) + y,,,d<H{R)a"-^m) - 2d^d'H{R)d,a”-^-'g(R)) ^/\^\ d"x 

^t^Ju 

+ f (R,v@'(R)n"'W(R)-R,v(@'(R)n'"H(R)))6/''^d”x. (29) 

Jm 

Using the equation (fT2ll we obtain the last integral in the above formula. Finally, one can put everything 
together and obtain 


6S1 =Ii + l2 + f’/J„ = -q f v'W(R)!F(n)@(R)6/’'^^d”x+ f {r^^O - 6/^ ^ d'b 

^ ^ Jm Jm 

, 00 n-1 „ 

+ X y /„ y (g,v (j“n'’H(R)Jan"-i-'@(R) + □'<K(R)n"-'@(R)) 


n=i ;=o 


- 2J^n'<K(R)Jva"“^“'@(R))5/’^ d”x. 


(30) 


□ 


Theorem 2.1. The variation of the functional lO is equal to zero iff 
~ 5 V'’H(R)!F(n)@(R) + (R.v'h - 

^ 00 n-1 

+ 2 E /« £ (^f,vJ“n''H(R)J,,n”-i-'@(R) - 2Jf,n'<H(R)Jva"“^“'@(R) + g,,vd'H(R)a"-‘g{R)) = 0. (31) 

n=l 1=0 

Proof. Since we have 5S = + 6 S 1 the theorem follows from Lemmas l2.3l and l2.6l □ 
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3. Signature (1,3) 

In the physics settings, where functional S represents an action, theorem 12.11 gives the equations of 
motion. From this point we assume that manifold M is the four-dimensional homogeneous and isotropic 
one with signature (1,3). Then the metric has the Friedmann-Lemaitre-Robertson-Walker (FLRW) form: 

+ r^de^ + sin^ Sd^^j. (32) 

Theorem 3.1. Suppose that manifold M has the FLRW metric. Then the expression dSTIl has two linearly independent 
equations: 

- 2<H{R)T{a)m) + {R<^ + 3n‘5) 

IotiG 

oo n-1 

+ + 2n'<H(R)n”-'@(R)) = 0, (33) 

«=i /=o 


ds2 = -dt2 + fl(t)2|_^ 


- \goo^{R)r{a)m) + (KooO - KooO) 

.. oo n-1 

-5E/.E (^oo<?“n''H(R)d„n”-i-'@(R) - 2doa‘'H{R)doa”-^-‘giR) + ^ooa''H(R)D”-'@(R)) = 0. (34) 

n=l 1=0 

Proof. The FLRW metric satisfies R,„, = jSft^v and scalar curvature R = ^ only on t, 

hence equations dSTIl for p v are trivially satisfied. On the other hand, equations with indices 11, 22 and 
33 can be rewritten as 

- mR)nn)g{R) + -0 + J^f„J^ (d“n'<H(R)d,n”-i-'@(R) + □'<K(R)n”-'@(R))) = 0. 

n=l 1=0 

Therefore these three equations are linearly dependent and there are only two linearly independent equa¬ 
tions. The most convenient choice is the trace and 00-equation. □ 

Corollary 3.1. For TT{R) = RF and @{R) = R‘^ the action ([2) becomes 

s = fjj^ + R’^mR^) ^ d”x, (35) 

and equations of motion are 

^(G,v + A^,v) - ^g,.R>’T{m‘^ + {R,v<^ - Krv^) 

1 “ (36) 

+ 2 £ Z (^f,vd"n'RPd,,n"-i-'R^ - 2d,,a’RPd,a^-^-‘R^ + g^,,dR’’n’’-'R^) = 0, 

n=l 1=0 

where O = pRP~^T'ia)Ri + qR‘>~^'F'(n)Ri’. 

Corollary 3.2. For 9L{R) = RF and @{Rf = R'^ the equations of motion (l36ll are equivalent to the following two 
equations: 

, oo n-1 

(4A R) 2R'’r(a)R‘i + (RO + 3nO) -i- E/«E (df'n'RPd^n"-i-'R'? + 2D'RPn"-'R'?) = 0, (37) 

n=l 1=0 
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167ig (*^oo ^9oo) - 29 ooR^'^{o)R‘^ + (Rqo® - Roo®) 

^ oo «-i (38) 

+ + goodRPa^-^R^ - 2dodKPdoa"-^-’R^) = 0 . 

n=l 1=0 


4. Constant scalar curvature 

Theorem 4.1. Let R = Rq = constant. Then, solution of equations of motion (l36ll has the form 


1. For Rq > 0, a{t) = -y ^ + ere ^ ^ + xe f 3 f lohere 9k^ - R^ot, cr, t e ]R. 

2. For Rq = 0, a{t) = yj-kfi + at + t, where + 4fcT = 0, 0 , t e ]R 


3. For Rq < 0, a{t) = \j ^ + o cos -J^t + t sin where 36k^ = R^(o^ + t^), a, t e IR, 


where k is curvature parameter. 
Proof. Since R = Rq one has 


/fl /a 

x2 k X 

'(—^ “ 

■ + ^ 

'fl 'fl 

' a^' 


;o. (39) 

The change of variable b{t) = a^(f) yields second order linear differential equation with constant coeffi¬ 
cients 


3h - Rob = -6k. ( 40 ) 

Depending on the sign of Rq we have the following solutions for b{t) 

6k [fof _ /Isf 

Rq > 0, b(t) - — + ae^ ^ ^ + Te ' ^ , 

Ro 

Rq = 0, b(t) - -kf + at + T, (41) 

6k /—Rq / —Rq 

Rq < 0 , b(t) = — -f- 0 cos y + T sin y 
Putting R = Rq = const into (l37t and ll38t one obtains the following two equations: 

+ = /oRr’(5Ro + (P + <()M 

Equations ll42ll will have a solution if and only if 

R'’^’"\Ro + 4 Roo)(Ro + ( 2 A - Ro)(p + q)) = 0 . ( 43 ) 

In the first case we take Rq + 4 Roo = 0 that yields the following conditions on the parameters 0 and t: 


Rq > 0, 9k^ = Rq0t, 

Rq = 0, 0^ + 4kT = 0, (44) 

Rq < 0, 361^ = Ro(0^ ■+ T^). 


□ 
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Solutions given in ||4T]| together with conditions ll44ll restrict the possibilities for the parameter k. 
Theorem 4.2. 1. If Rq >0 then for k = 0 there is solution with constant Hubble parameter, for fc = +1 the 

solution is a{t) = .^^cosh j ^ and for k = -1 it is a{t) = sinh j ^ , where 

cT + T = -^ cosh (p and a - t = ^ sinh (p. 

2. IfRo = 0 then for k = 0 the solution is a(t) = yfr = const and for k = -1 the solution is aft) = |f + ||. 

, where o = ^ cos <p and 


3. If Rq < 0 then for k = -1 the solution is a(t) = 
^sincp. 


cos 


2 ( - <P) 


Proof. Let Rq > 0. Set k = 0 then we obtain solution with constant Hubble parameter. Alternatively, if we 
set A: = +1 then there is cp such that o + t = ^ cosh cp and o - t - ^ sinh (p. Moreover, we obtain 


w ^ 12 , 2 1 /Ko 

6(t) = -^cosh 2 


12 1 
a{t) =A/^cosh- 


At the end, if we set A: = -1 one can transform b(t) to 


12 . ,2l 

b{t) = — smh^ - 
i<0 z 


t (p 


a{t) = 


sinh ■ 




t + cp 


(45) 


(46) 


Let Rq = 0. If A = 0 fhen fimcfion b{t) and consequenfly a{t) become consfants which leads to a solution 
a{t) = xfr = const. On the other hand if A 0 fhen we can wrife 


b{t) = ~k{t--^)f 


(47) 


If A = +1 then there are no solutions for the scale factor a{t), because b{t) < 0. On the other hand, when 
A = -1 the scalar factor becomes 


fl(f) = |f+-|. 


(48) 


In the last case, let Rq < 0. If A = -1 we can find cp such thaf o = cos cp and t sin^ and rewrife 


a(t) and b(t) as 


-12 2 I / Ro, 




1 

COS - 
2 


Ro, 


(49) 


In the case A = +1 one can transform b(t) to b(t) - ^ sin^ j 
yields no solutions. □ 

Theorem 4.3. If in ll43| we take 


-■y f - (pj, which is non positive and hence 


r: 


p+(|-i 


(Ro + (p + (?)(2A-Ro)) = 0 


(50) 


then: 


1. For p + q>l there is obvious solution Rq = 0. In particular ifp + q = l then ll50t is satisfied for any Rq 0 if 
A = 0. 

2. For p + q = 0 there is no solution. 

2A(p + q) 

3. For p F q + 9,1 there is a unique value Rq - — - -^ that gives a solution. Since p and q are integers the 

value of Rq in the last equation is always positive, and for k = 0 the solution b(t) is a linear combination of 
exponential functions. 

Proof of fhis theorem is evident. 
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